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Abstract. 

The dynamics of quantum systems strongly depends on the local structure of 
the Hamiltonian. For short-range interacting systems, the well-known Lieb-Robinson 
bound defines the effective light cone with an exponentially small error with respect to 
the spatial distance, whereas we can obtain only polynomially small error for distance 
in long-range interacting systems. In this paper, we derive a qualitatively new bound 
for quantum dynamics by considering how many spins can correlate with each other 
after time evolution. Our bound characterizes the number of spins which support the 
many-body entanglement with exponentially small error and is valid for large class of 
Hamiltonians including long-range interacting systems. To demonstrate the advantage 
of our approach in quantum many-body systems, we apply our bound to prove several 
fundamental properties which have not be derived from the Lieb-Robinson bound. 
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1. Introduction 

The fundamental features of quantum many-body systems are strongly restricted by 
the local nature of Hamiltonian. Such restrictions give us a lot of useful information 
in analyzing universal properties of matters. One of the prominent examples is the 
Lieb-Robinson bound [1, 2], which characterizes the velocity of information propagation 
in non-relativistic quantum systems; in other words, we can define an approximate 
“light cone” with an exponentially small error. Based on the Lieb-Robinson bound, 
we can grasp fundamental restrictions to quantum dynamics: entropy production rate 
after quench [3, 4, 5, 6], entanglement growth [7, 8, 9, 10, 11, 12, 13, 14], complexity 
of quantum simulation [15, 16, 17, 18, 19], and so on. Moreover, the Lieb-Robinson 
bound also provides us powerful analytical tools to give foundations of quantum 
many-body systems, from condensed matter physics to statistical mechanics: Lieb- 
Schultz-Mattis theorem [20, 21], exponential decay of bi-partite correlation [22, 23, 24], 
stability of topological order to perturbation [25, 26, 2, 27, 28], quantization of the 
Hall conductance [29, 30], thermalization problem [31, 32], equivalence of the statistical 
ensembles [33], etc. In these results, the locality of interactions plays essential roles. 
Thus, the principle of locality has shed new light on our understanding of fundamental 
many-body physics. 

More recently, with the progress of experimental technology [34, 35, 36, 37], there 
has been considerable interest in the potential of the locality analysis in long-range 
interacting systems, both from theoretical [38, 39, 40, 9, 41, 42, 43, 44, 45, 46, 47] and 
experimental [36, 37] viewpoints. In such systems, we can also define the approximate 
light cone as in the case of the short-range interacting systems. However, the light- 
cone is usually nonlinear to the time except some special cases [42], and moreover the 
transport of information can be bounded only polynomially [41, 46, 23, 36, 42] with 
respect to the spatial distance outside the light cone. The primary reason is that the 
Lieb-Robinson bound focuses on the velocity of the information transfer, whereas the 
long-range interacting systems can transport information immediately in principle. In 
this way, the causality allows us to analyze the system in the looser way in comparison 
with the case of the short-range interacting systems. This indicates that we may not 
grasp all the restrictions due to the locality of the Hamiltonian in terms of only the 
spatial distance. Here, we use the term of “locality” in more broader meanings [40] 
as distinguished from the spatial locality, i.e., to what extent can quantum systems be 
described by a collection of local degrees of freedom, which are only loosely correlated 
with each other? 

In the present paper, we give a qualitatively new bound for dynamical properties 
in terms of the number of spins instead of the spin-spin distance ; roughly speaking, 
we focus on how many spins can correlate with each other after time evolution. In 
order to make our concept clear, we first consider a classical process in which a source 
transfers information to receivers. We now assume that each of the elements sends 
information to any other 2 elements per unit of time (see Fig. 1). Then, the number of 
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Figure 1 . Schematic picture of classical information spreading. We consider a 
classical system with N elements and assume that each of the elements is connected 
to arbitrary ones. We start from one information source (e.g., coded by 0 or 1) and 
consider the process that each element which has received the information can send 
the information to other 2 elements per unit of time. For example, after a unit of 
time, 2 elements receive the information and then totally 3 elements becomes the next 
information source. In the same way, after two units of time, totally 9 elements share 
the information. Repeatedly, the number of receivers after n units of time, say N n , 
increases as 3". By analogy, we also expect that the information spreading of spins 
should be similarly suppressed in the quantum many-body dynamics. 


o 

t = o 
(N 0 = 1 ) 


elements which share the information can be bounded from above by e 0 ^\ Our main 
purpose is to give a quantum version of the bound in the form of an operator inequality. 
As we shall see shortly, our new bound characterizes the number of particles which 
support many-body entanglement such as the topological order [2] and the macroscopic 
entanglement [48, 49]; this indicates that the global entanglement induced by a small¬ 
time evolution should be suppressed up to exponentially small error. It allows us to 
obtain novel strong properties on various kinds of fundamental physics in quantum 
many-body systems which have not be derived from the Lieb-Robinson bound. 

2. Model and formalism 

We consider a spin system of finite volume with each spin having a d-dimensional 
Hilbert space and label each spin by i = 1, 2,... N. We denote the set of all spins 
by A := {1,2,..., A/”} We denote partial sets of sites by X, Y, Z and so on and 
the cardinality of A", that is, the number of sites contained in A", by |A| (e.g. 
X = {ii,i 2 ,..., i|x|})- We here define the q-local operator as follows: 

r(9) = (g-locality), (1) 

X\<q 

where each of the {yx} is supported in a finite set X C A. In other words, the q- local 
operator contains up to g-body coupling. 

Here, we assume systems which are governed by fc-local Hamiltonians with k = 
0 ( 1 ): 

H=Y, h A'- 

|X|<k 


( 2 ) 


4 


We assume the time-independence of the Hamiltonian for the simplicity, but the 
discussion can be generalized to the time-dependent Hamiltonian H(t). Note that we 
make no assumption on the geometry of the system, and the coupling can be arbitrarily 
long ranged. Instead, as a normalization factor, we introduce the parameter g of the 
Hamiltonian: 

11 hx 11 < g for Vi G A (g-extensiveness) (3) 

X:Xli 

with || • ■ ■ || the operator norm (i.e., the maximum singular value of the operator); we 
refer to that the Hamiltonian is ^-extensive if it satisfies the condition (3). This implies 
that the energy associated with one spin is bounded by a finite value g. Note that the 
norm of the Hamiltonian ||ih|| increases at most linearly with the system size N, namely 

N N 

<X>*II^£ E UM< = (4) 

X i= 1 X:X3i i= 1 

We notice that the class of h-local Hamiltonians covers almost all realistic quantum 
many-body systems not only with short-range interactions but also with long-range 
interactions. 

For the basic analysis of the fc-local Hamiltonian, we utilize the following theorem: 



Theorem 2.1 Let H be a k-local g-extensive Hamiltonian and be a q-local operator; 
note that the T^ may not be extensive as in ( 3 ). Then, for an arbitrary positive integer 
q, we can obtain 

||[tf,r (9) ]|| < a|||T (,?) || with A:=6 gk 2 . (5) 

Note that the operator [H, T^] is still at most (k + g)-local. In the case where T^ is 
supported in a local subset Z C A, namely = j z (\Z\ — q), we can simply prove 
the theorem as follows: 

II [H,3z \|| < \\l h xi3z}\\ < 2 W hx W ' \M\ ^ 2 9\ z \ ■ hz\\ < 6^gfc||7z||-(6) 

X-.XnZj4 i&Z X:X3i 

When rw is a general g-local operator, however, the proof cannot be given in a simple 
way but a bit technical. To show the point, we expand a g-local operator T^ as 

r (,) = ;r 7 z. 

\Z\<q 

Now, the difficulty lies in the fact that we cannot utilize the following simple estimation; 
we have 

l|[ff,r«]||< ^ \\[H, lz \\\<2 hzll, 

\Z\<q \Z\<q 

whereas we cannot generally ensure 

E Itell«lir (,) ||. 

\Z\<q 
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For example, let us consider the following 2-local operator 

r<2) = ^E^.*®4 

i<j 

where are uniform random numbers from —1 to 1. For this operator, we can 

obtain Y2i<j II ® a j W/N = O(N), but ||r (2) || = 0(VN), and hence 

ll[-ff.r <2| ]|| < \W H M\\ < ll7ij y ^ = 4 9 0(Viv)||r< 2 >||, 

\Z\<2 i<j 

where the second inequality comes from (6). This is much looser than the inequality (5) 
which conies from Theorem 2.1. We give the full proof in Appendix A. 


3. Main results 


In order to mathematically apply the classical discussion on the information sharing 
to quantum cases, we consider the dynamics of the /c-locality of operators. We 
initially consider a go-local operator and investigate its time evolution: T^ qo \t) = 
e -tHt-p(.qo) e iHt _ After a time evolution, the operator T^(t) will be no longer a g 0 -local 
operator but may be approximated by another g-local operator with q > g 0 , say 
We now regard go an d g as the numbers of particles which share the information at the 
times 0 and t as in Fig. 1, respectively. We then expect that the approximation can be 
rapidly improved beyond q > q 0 e°^ from the classical discussion. Indeed, we prove 
the following theorem for the minimal error of the approximation: 


Theorem 3.1 LetU(q) be a set of q-local operators and consider an arbitrary q 0 -local 
operator 6 lA(qf). Then, for an arbitrary real-time evolution with t > 0, there 
exists a q-local operator which approximates the operator (t) with the following 
error: 


inf (||r (9o) (f) -rj 9) ||) < 8||r (9o) ||rft;t] 


r \ q) eu{q) 


exp 


1(1 


% 


( 7 ) 


with k = 4A, £ = k/ log 2 and r t = 2 ^ _ where [• • ■] denotes the ceiling function. 
The same inequality holds for t < 0 by replacing t with |t|. 


Because the function r t increases as e 0< ^\ the time-evolution of r( 90 ) can be well 
approximated by a (goe°^))-local operator. The upper bound is meaningful as long as 
t < logA^, which is qualitatively consistent with the threshold time of the breakdown 
of the Lieb-Robinson bound for long-range interacting Hamiltonians. As in the case of 
the Lieb-Robinson bound [41, 42], we might improve the present theorem by explicitly 
considering a spatial structure of the system, for example, the power-law decay of 
interaction. 

We also mention the case where the system is governed by a short-range interacting 
Hamiltonian. In this case, we can obtain much stronger restrictions [2], Let us consider 
an operator T L which is supported in a region L. After a short time, the operator F L (t) 
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Figure 2. Schematic picture of the proof of Theorem 3.1. For small-time evolutions, 
we can obtain the approximation of T^ qo \t) directly by the use of the time expansion 
as in Eq. (9). For long-time evolutions, we split the total time range [0, t] into n 
intervals with a length St. In each of the intervals, we can apply the result for small¬ 
time evolutions, and hence we can connect the approximations as follows; we first 
approximate T^ q °\St) by second approximate r^^((5t) by T^ and so on. Then, 

we can finally prove that the operator satisfies the inequality (7) by choosing the 
set {g m ,r^ m) }™ =1 appropriately. 


is no longer supported in the region L, but is approximately supported in some region 
having distance l from L. The Lieb-Robinson bound ensure that the accuracy of this 
approximation becomes precise exponentially as the distance l increases beyond 0(t). 
In other words, the support of T [ft) enlarges at most as Poly(t) instead of e°^\ 

Proof of Theorem 3.1. For the proof, we first obtain the upper bound for a small-time 
evolution: 

mf air<®>(() - r<»>io < 2«/‘. ( -f )( :;; ,/ V ( - ) ii m 

r ( t q) eu( q ) l-Kt/2 

for t < 2 /k. In the derivation of (8), we use the direct expansion of e~ lHt T ( ' qa ' ) e lHt 
according to the Hadamard lemma of the form 

m 

r (w) (f) = E [H ' [H ' "' {H ’ r(so) ]]" '1 ■• (9) 

m =0 

We terminate the expansion (9) at m — mo so that the expanded operator may be 
g-local. We then estimate the error due to the termination to prove the bound (8). We 
show the full proof in Appendix B. 

We, however, cannot utilize the expansion (9) in order to obtain a meaningful bound 
for t > 2 /k. In obtaining the inequality (7), we will have to utilize the fact that e~ lHt is 
unitary f For this purpose, we split the time range [0, t) into n intervals (Fig. 2) such 
that t/n < 1 /k. We here denote the length of the interval t/n by St: 

St := — < — with n := Int I. (10) 

n k 

We also define t m := mSt for m = 0,1, 2, ... n with t n = t. Note that in each interval 
we can now apply the upper bound (8). 

f Note that the bound (8) can be also applied to the imaginary time evolution without 

the unitarity condition. The approximation by the finite expansion of e Ht T ( - qo ^e~ Ht will becomes less 
accurate beyond a certain time t c , which comes from the fact that the norm of e Ht T^ q °^e~ Ht rapidly 
increases for t > t c . Without the unitarity, we cannot arrive at the inequality (7) from (8). 
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In the following, we connect the approximations of T ( qo >(t) from the first interval 
[0, ti) to the last interval [t n _ 1 : t n ]. We first approximate the time evolution T^^St) 
with a gi-local operator r^. We second approximate T^(<5f) with a g 2 -local operator 
By sequentially repeating this process, we define a set of operators {T 
so that they may approximately satisfy ||r^^(5£) — ~ 0, respectively, where 

q n < <?• We thus obtain the approximation of r^°^(t) by the use of the set {r[®'™' l }” =1 : 

l|r<«" ) - r (®)(()|| 

n 

= lEibrbn - u - ifc ,l (i. ~t m +«)]| 

m= 1 
n 

< EiihrA - t m ) - - t m +*>ii 

m= 1 
n 

= Eiibr’-brEwn, in) 

m= 1 

where we used the equality t m+ \ — t m — St in the first line, the norm invariance during 
the time evolution in the third line, and we set =: T^b Note that we use the 
unitarity of the time evolution in the third equality. 

We can then show the appropriate choice of the set {q m , Prom the 

inequality (8), we can prove that there exists a set {q m , T^ m ' ) }^ =1 such that: 

l|rtr } - rfc l} ( 5 £)|| < A(A + i) m - 1 ||r (9o) || (12) 

for m — 1,2,, 7i, respectively, with A := 4exp[—| ^ — g 0 ^]. The derivation is given 

in Appendix C. 

By combining the inequalities (11) and (12), we have 

n 

||rg n) - r {qo \t)\\ < A ( A + i) m_ 1 l|r (9o) || 

m =1 

= [(A + i) n -i]-||r (,?o) ||. (13) 

We can always find an operator T such that ||T^ n ' ) — r^°^(£)|| < ||r^ 0 )|| (e.g. 
r)'( n = 0), and hence we only have to consider the range (A + l) n — 1 < 1 in the 
above inequality and obtain 

(A + l) n - 1 < 2 1-lA • n • A < 2nA. (14) 

The inequality (13) reduces to the inequality (7) due to the inequality (14) and the 
definitions of A and n. This completes the proof of Theorem 3.1. □ 

§ Proof of (14): the inequality (x + l) n — 1 < 1 is satisfied for x < 2 1 /" — 1. By the use of the fact that 
(x + l) n — 1 is the concave function for x > 0, we have {x + 1)" — 1 < nx(xo + l) n_1 for 0 < x < xg 
with xq a positive constant. By choosing xq = 2 1 /™ — 1, we obtain (x + l) n — 1 < nx 2i" _1 h n . □ 




4. Several implications 


4-1. Stability of the topological order 


We here prove that the topological order is stable after the time evolution over t < log N. 
For the definition of the topological order, we follow the same discussion as in Ref [2], 
The concept of the topological order is usually defined with respect to Hamiltonians 
rather than quantum states. However, we have several common properties which the 
topological ordered phases always satisfy. 

Slightly generalizing the definition in Ref. [2], we here define that a quantum state 
| if) exhibits the topological order if and only if there exists another quantum state \if) 
which satisfies 


(V>|r (q) | fj) = <' 0 |r (9) |^> and (V’|r (9) |^) = o 


for arbitrary g-local operators with q = 0(N q ) and q > 0. Here, the main difference 
from Ref. [2] is that we apply generic g-local perturbation instead of the spatially local 
perturbation. If the quantum state satisfies the property, the coherence between the two 
states \if) and \if) can never be broken by any kinds of local operators. It is known that 
topologically ordered phases usually satisfy these conditions, for example the ground 
states of the Kitaev’s toric code model [50]. I 11 Ref. [2], for short-range interacting 
systems, they proved by using the Lieb-Robinson bound that the topological order 
persists at least for t < iV 1//£> ( D: system dimension). This is contrast to the present 
case where we apply the new bound (7) and ensure the time of the stability for t < log N. 

In considering the stability of the topological order, we define the topological order 
with error (g, e q ) as 


\{fj\T (q) \fj) - (^|T (9) |^)| < e q and \{^\T^)\ < e q , 

where ||T^|| = q. For exactly topologically ordered states, we have e q = 0 for 
q < g 0 = 0(N a ) with a > 0. We will see that after small-time evolution the error e q 
can be small sub-exponentially with respect to the system size. In evaluating the error 
(g, e q ) in the case of the fc-local Hamiltonians, we apply the bound (7) in Theorem 3.1 
instead of the Lieb-Robinson bound. From the similar discussions as in Ref. [2] , we can 
obtain 


e q {t) < 2 \nt] exp 



J 


where we use Theorem 3.1 in evaluating |(^(t)|r^|^(t)) — ('0(f)|r^|'0(f))j = 

|(V>|r^(— f)|V>) — ( , 0|r^(— £)|-0)|. If we consider q < q 0 /[2r t ], the error e q is bounded 
from above by 0{t ) - exp [— qoe~°^]. This means that the error of the topological order 
is small exponentially with respect to go as long as t < logg 0 = O(logA^). 


4-2. Dynamics of probability distribution of macroscopic observables 

We finally discuss how to connect our Theorem 3.1 to observable quantities. As an 
example, we consider an upper bound for distribution function for extensive quantities 
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A such as A = JA =1 a* with 1111 = 1 for i e A. Let us consider a product state |Prod) 
and its time evolution e~ lHt \Prod). Initially, the spins are independent of each other and 
the distribution of A for |Prod) is concentrated with the standard deviation at most of 
0(y/N) due to the Chernoff bound [51]. After the time evolution, the spins can couple 
with each other but still maintains the local independence approximately. Hence, we 
expect that the distribution of the operator A should be still concentrated in |Prod(£)). 
Indeed, we can prove the following theorem: 

Theorem 4.1 Let Tif z be a projection operator onto the subspace of the eigenvalues of A 
which are in [z, oo). Then, the spectrum of A in |Prod(£)) is exponentially concentrated 
as 

ll n >(A>+/?l Prod W)ll < c i ' exp(— - J 1 — ), (15) 

with Ci and c 2 0( 1) constants, where ( A) is average value with respect to |Prod(i)). 

After a short time t = 0(1), the distribution of A is still strongly concentrated 
with a standard deviation e°^\f~N. Thus, in the state |Prod(£)), spins are still locally 
independent of each other. We note that this also implies no macroscopic entanglement 
in terms of the quantum Fisher information [48, 49]. We expect that the inequality (15) 
can be improved to the Gaussian form by using the recent technique [52], 

Proof of Theorem f.l. For the proof, we focus on the fact that the product state |Prod) 
is given by ground state of a 1-local Hamiltonian, say H p : H p = YliLi hi with ||/q|| = 1 
for i G A and // p |Prod) = —A|Prod). Note that the spectral gap A E between the 
ground state and the first excited states for H p is 0(1). Then, the state |Prod(i)) is 
also a gapped ground state of H p (t) = hi(t). We now expand H p (t) by the use of 

n A •= n A - 

[xrt,rt+xrt) ' x m 

H P (t) = E nin P mi 

x,x' 

■=^^^(x,s\H p (t)\x',s')\x,s)(x',s'\, (16) 

x,x' s=l s'=l 

where we denote n A = ^Js=i l x > s )( x > s l with d x = rank(n A ). 

By applying Theorem 3.1, we obtain 

||n^p(t)n^|| < C v e-^ x - x '\ (17) 

with 

C v = 8Ne 5/ W\Kt], n = 1/(20, ( 18 ) 

where the derivation is given in Appendix C.l. This way, we can formally regard 
the Hamiltonian (16) as a tight-binding Hamiltonian; the position x corresponds to 
the eigenvalue of A. Remembering that the state |Prod(f)) is the gapped ground 
state of H p (t), the distribution of the position x should be localized due to the 
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spectral gap A E [53, 54, 55]; from Ref [55], the localization length is proportional 
to ^C v/ i 3 /AE oc yftN with A E = Oil). Now, 1 jr t times of the eigenvalue of A 
corresponds to the position x, which yields the localization length of A which is smaller 
than 0(rtVtN). This complete the proof. □ 

Before closing this subsection, we discuss the relationship to the spin squeezing [56] , 
where the magnetization along a certain axis can be squeezed by broadening the variance 
along another axis. The squeezing protocols contains the process to create large 
fluctuation and hence Theorem 4.1 is applicable to the estimation of necessary time 
for the squeezing creation; for example, let M x , M yi and M z be the usual collective spin 
operators and assume that the total spin amplitude is O(N). Then, the uncertainty 
principle ensures AM x AM y > CN with C a constant of Oil) and A(-) denoting 
the fluctuation. Theorem 4.1 implies that the necessary time for the squeezing of 
AM X /AM y = £ sq < 1 (i.e., A M y > y/CN/£ sq ) is at least proportional to log(^ 1 ). 
This time scale is apparently not consistent to the many previous works [56], in which 
the necessary time is Oil) even if £ sq = N~ a with a > 0 for protocols using two-body all- 
to-all interactions like Lipkin-Meshkov-Glick model [57]. The key point is that we now 
assume the extensiveness of the Hamiltonian (3), whereas in spin squeezing literatures, 
the Hamiltonians are super extensive. By taking this point into account, we can resolve 
the inconsistency. 

5. Outlook 

We have given a new bound on the quantum dynamics which are governed by the fc-local 
Hamiltonian and shown some applications which have not be derived from the Lieb- 
Robinson bound for long-range interacting systems. Our main theorem characterizes 
the number of spins which cause many-body quantum effect due to time evolution. 

As further applications, it is one of the most important problems whether we can 
apply the present results to the quasi-adiabatic continuation [58, 27] for the /c-local 
Hamiltonians. In more detail, let the Hamiltonian be H{s) — H + sV with 0 < s < 1, 
where H and V are fc-local Hamiltonians, respectively. Here we assume that the ground 
state of H(s), say |i?o(s)), is non-degenerate and is separated from the excited states 
by a non-vanishing gap for 0 < s < 1. Then, the evolution of the state can be described 
similarly to the time evolution: f^\E 0 (s)) = iD(s)\E 0 (s)), where D(s) is defined as 
the adiabatic continuation operator. We therefore can conclude that the parameter 
evolution of lE'o(s)) is formally equivalent to the time evolution by D(s). Our interest is 
whether the operator D(s) can be approximated by fc-local operators. In that case, we 
can analyze the perturbative effect to the ground state by the use of Theorem 3.1, from 
which we can generalize the stability analysis for the topological order in the short-range 
interacting systems [25, 26, 28]. 

Another interesting direction is to identify an efficient description of the state 
|Prod(t)), for example, by the use of the tensor network state [59]. Because the 
multipartite effect is highly suppressed as long as t = 0( 1), we expect that the small 
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time-evolution can be efficiently simulated. Moreover, if the Hamiltonian contains 
randomness in its interaction, the locality of the multipartite coupling might persistently 
maintain beyond the time scale t > log IV, in the similar manner as the short-range 
interacting systems [7], where the entanglement growth is logarithmically slow. 

Finally, can we observe our bound experimentally? In particular, Ref. [36] has 
demonstrated the Lieb-Robinson bound in long-range interacting systems. Our new 
bound in Theorem 3.1 can be in principle observed in the same experimental setup, for 
example, by looking at the probability distribution which follows Theorem 4.1. 
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Appendix A. Proof of Theorem 2.1 

We show the proof of Theorem 2.1. For the proof, we take the following two steps. 
(Step 1) We first consider a class of the commuting Hamiltonians H c as 



(A.l) 


|X|<fc 

and prove the inequality 


llpf c ,r ( »>]|| <6g ? l|r ( »>||. 


(A.2) 


(Step 2) Secondly, we prove that any /.-local ^-extensive Hamiltonian can be 
decomposed into the sum of commuting Hamiltonians: 



(A.3) 


m=l 

where each of {H^} is /t-local and (gfc)-extensive. 
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After these two steps, we can prove Theorem 2.1 as 

1 ” 

U 171= 1 

1 " 

< F ^6^g||r (,) ||=6^g||r (5) ||. (A.4) 

m= 1 

This proves Theorem 2.1. In the following subsections, we will prove the statements in 
the Steps 1 and 2. 

Appendix A. 1. Step 1 

We here prove the upper bound (A.2) for || [ H c , T^] ||. We first decompose H c as follows: 
H c = H' c + 6H C , 

where 

OO 

H' c := ^ e(j + l/2)U Hej+e) , 5H C := H c - H' c 

j=—o o 

and II| {J)£J+£ ) is a projection operator onto the eigenspace of H c with the eigenvalues 
[ej, ej + e). We set the value of e afterward. Note that the operator II[ e j )( y +e ) may be the 
null operator. From the definition, we have 

ll^ll < 

We then obtain 

||[T^, H c ]\\ = \\[T^,H' C + 6H C }\\ < ||[T (9) ,i/ ,c ]|| + ||[T^, 5H C ]\\, (A.5) 

which necessitates that we calculate ||[T^, H' c ]\\ and || [T^\ 6H C ] || separately. 

We first obtain the norm of [T^ q \SH c ] as follows: 

||[r {9) ,«T]|| < 2||T^|| • \\SH C \\ < e||T^||. (A.6) 

We second obtain 

3,3' 

= ™ i)n[ e yej+e)r^n[ e j/ )e j/ +e ). (A.7) 

3,3' 

Because we can obtain the norm of [T^, H' c ] from the equality 

||[ry i/ c ]|| = max |(^|[r*y H c \\^)\, 
m 

we, in the following, calculate the upper bound of |(^|[T^, H' c ]\^)\ for arbitrary 
quantum states |^). From Eq. (A.7), we have 

-j)Wi n M«+.) r<s)n i«'«'+.)i^) 

3 , 3 ' 
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— 'y ' e \j j I ■ ||<^|n [ejW+e) || • \\u [ejtej+e) r^u [ejl>ejl+e) \\ ■ \\u [ejl}ejl+e) 






3,3 


where aq := \\U [ejtej+e )\ip)\\ and := ||n [e j iei+e )r^II[ e j/ |ej v +e) ||; note that J2j a j = 1 - 
Because we are now considering the commuting Hamiltonian, due to g- 
extensiveness (3), we have [38, 40] 

]|n [£ ;/, 0 o) r { 9 ) n ( _ 0OijB] || < ||r^|| for \e' -e\< 2 gq , 

||n [S /,oo)r ( ' ?) n ( _ 00iE] || = 0 for \E'-E\>2gq, (A.8) 

where the first inequality comes from the trivial bound of Un^o^T^ n(_oo ,£]ll < 
lin^.oo) || • ||r (,?) || • i|n ( _ 0O)B ]|| = ||r (y) ||. This gives the following inequality: 


rfj, < ||r<«>|| for \j'-j\<i + 2 f-, 


T% = 0 for + 


2 gq 


3,3 


Because of the inequality (A.9), we have 

|(#rw,= £ V 


|j'-j|<l+29<j/ e 


<e||r (,) ll £ U'-J 

<£||r (,) ll £ 

\j'-j\<l+2gq/e 

[1+2|2J 

< e||r (,) |l £ \f\ J2 

i'=- LI+^J 


4 + a 3 


O'■ 


< e||r (,) ll 


1 + 


2^g 


e J 


1 + 




e J 


+ 1 


where [■ • -J denotes the floor function. 

The inequalities (A.6) and (A. 10) reduce the inequality (A.5) to 


, 2 gq 
1 H—— 

e J v l 


2 gq 
1 H —— 
e J 


+ 1 


||[r(’>,f/ c ]||< £ ||r<’)||+ £ ||r'«'|| 

We here choose e = 2 gq + 5e (<5e > 0) and obtain 
||[r ( 9 o) ,iT]|| < (6gq + 35e)||r^°^||. 

By taking the limit of 5e —> +0, we finally obtain ||[T^,i7 c ]|| < 6^g||T^||. 


(A.9) 


(A.10) 


Appendix A.2. Step 2 

We now prove the existence of the decomposition (A. 3). For the proof, we first introduce 
a parameter e and define unit operators {hx} as follows: 

h x 


h x ■= e 


II h 


x 


(A.ll) 
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where {hx} are components of the Hamiltonian. By the use of { hx }, we can dehne the 
‘discretized’ Hamiltonian H' as 


H' = Y, N hx~hx 

X 

I H' 


with 


N hx := 


II h 


x\ 


(A. 12) 


Note that we have \\H' — H\\ = 0{eN ), which vanishes in the limit of e —>■ 0. Because 
of the extensiveness of the Hamiltonian, the number of unit operators {hx} which one 
spin can contain should be bounded from above by 


(A. 13) 

X3i 

We, in the following, decompose the Hamiltonian H' into n commuting 
Hamiltonians by the use of {Ay}, where n is an integer. We here correspond 

a set of {h (i)}^ i to one commuting Hamiltonians ifb, namely 

-X-m 


H} n ~ h Y (i) + h ( 2 ) + • • • + h y (N m ), 

m x\-m -A-m 

s.t. X$nXjQ = 0 for Vj,k (A. 14) 


Note that the Hamiltonian H c is commuting, fc-local and e-extensive because of 
\\h „(i) II = e. Hence, if we can decompose the total Hamiltonian H' with n = k\g/e I, we 

-A-m 

obtain 

1 

H‘ = r £ nH (AA5) 

m= 1 

The commuting Hamiltonians are fc-local and (gfc)-extensive. Thus, by 

taking e —> 0, we have \\H' — H\\ — y 0 and this completes the proof of the 
decomposition (A.3). 

In the following, we will prove that n = k\_g / e\ is a sufficient number of commuting 
Hamiltonians {H^} to construct H'. According to the definition of in Eq. (A. 14), 
we denote the support of H{ n by L m , namely 

L m = XW u A'if IJ ■ ■ ■ U Xg-K (A. 16 ) 

In the construction of {H^}, we decompose the Hamiltonian H' such that \L l \ < \L 2 \ < 
• • • < |An|- 

We hrst collect the units {hx} for so that \Li\ may be as large as possible, 
mathematically, 

H' Lt := £ N h x hx= 0. (A.17) 

x-.xeLf 

This means that there are no unit Hamiltonians {hx} outside of L\. Second, we also 
collect {hx} into so that \L 2 \ may be as large as possible: 

H = 0 with H {2) ■= H' - H[. (A.18) 

By repeating this process, we construct the commuting Hamiltonians {H^} so that they 
may satisfy 

m— 1 

HW = 0 with H^ m) := H' - V H c . 

/ ■> 3 

3 =1 


(A.19) 
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It means that there are no unit Hamiltonians {hx} outside of L m in considering 
Note that if L n+1 = 0 we also have = H' — ^” = L H/ — 0 and the complete 

decomposition has been achieved. In the following, we have to prove L n +i = 0 with 
n = k[g/ ej. 

For the proof, we assume L n+1 ^ 0, or equivalently H^ +1 ^ 0, and prove the 
contradiction. We use the fact that any terms {h v u) j+i" 1 supported in L n +1 should 
satisfy 

Xj£h n L m ^ 0 (A.20) 

for rn — 1,2,3,..., n, which comes from the condition (A. 19). Because of (A.20) and 
A^, < k, at least (n/k) subspaces in {L m }^ =1 have a common support; for example, 
if X^, = 1 or X^| contains only one spin (e.g. spin i ), the relation (A.20) means 
that all of {L to }^ =1 contain the spin i. 

Therefore, there exists a set such that 

hi L mi n L m , 2 n L m3 n • • • n L mn ^ k ^ 0. (A.21) 

We denote this common support by L. Then, the spins in L should be contained in all 
the Hamiltonians an d Hfi+i, whereas, due to the inequality (A. 13), one spin 

can contains up to (n/k) unit operators {hx}', note that [g/e\ = n/k. Thus, we prove 
the contradiction. □ 


Appendix B. Bound ( 8 ) for a small-time evolution 


We here prove the inequality (8) for a small-time evolution: in order to obtain the bound, 
we expand e ~ lHt Y ( - cl0 '>e lHt by the Hadamard lemma as in Eq. (9), which we reproduce 
here: 


e -itfip(go ) e iHt _ 


-ay 


771=0 


m\ 




(B.l) 


where g 0 = T^ 90 ) and L m = [H, [H, ■ ■ ■ [H, T^®)]] • ■ •]. This expansion can be terminated 
if the expansion converges rapidly as m increases. The termination at m = m 0 gives 
the local approximation of the operator r(®+), namely 


mo 


it y 


r (9) = E —r L - 

ml 


(B-2) 


m=0 


Because L m is at most (go + fcmo)-local, in order to make the operator (B.2) less than 
or equal to g-local, we take 

Q-Qo 


mo = 


k 


Our purpose now is to calculate the error due to the cutoff of the above expansion, 
namely 


E 

777=7770 + 1 


(-ay 


ml 


t r 

< E -riiM. 

* m I 


777=7770 + 1 


ml 


(B.3) 
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Now, we can apply Theorem 2.1 to evaluate \\L 
ill L 


m 


= || [H, [if,-••[#,['<*>]]■■■] || 

q 0 qo + kq 0 + 2k q 0 + (m~l)k 


^ k k k 

= A m (n +^- ! 1 ) ! ||r ( g ° ) || < 2 r+m m!A m ||r (9o) ||, 

where we define r := q 0 /k and [>] +m- 1 C m < 2^1 +m_1 < 2 r+m . We thus obtain 


k 


l|r (,s) || 


00 frn 00 ( c 2\t^ q - qo ' > / k 

—\\ L m\\<^ r (2At) m ||r (go) || < 2 q/k { ^_ ow ||r (go) ||, 


2\t 


(B.4) 


m=mo+l m=m o+l 

for t < (2A) _1 , where we used the inequality m 0 + 1 > (q — qo)/k in the last inequality. 
By replacing 2A = k/2, we can obtain the inequality (8). 


Appendix C. The proof of the inequality (12) 

We here prove the existence of the set {q m , T^ m ' ) }^ =1 which satisfies the inequality (12). 
Because of 5t < l//c, we can apply the bound (8) to estimate the norm ||T^ m ' ) — 
r fc l) (5t)||. As shown in the next subsection, we can find a set {q m , r^ T ^}™ =1 such 
that 

llhr ) -rSr_-; l) (*)ll<lirfc l) l|A 

for m = 1, 2,..., n, respectively, where 

1 f q 


(C.l) 


A := 4 exp 


f \n 


- qo 


(C.2) 


By the use of the inequality (C.l), we can obtain the following inequality: 

||r^ m) || < (A +i) m ||r (9o) ||. (c.3) 

We prove this inequality by the induction method. For m = 1, we have 

lir^H = \\r\f - rg o) + rg o) || < ||rf } - rg o) || + ||rJ 0 9o) || < (A + i)||r^||, 

where the last inequality came from (( !) and we used the definition of = T^°b 
We then assume the inequality (( A3) for m < m o and prove it for m = m 0 +1 as follows: 


I j-\(<7mQ + l) | 
I tmQ + 1 


{Quiq + i) -p(<7 mo) | -pi(^ m Q ) | 


I ^mQ + 1 


+ r 


< (A + i)||r£ no) || < (A + i) m o+i|| r (9i 


!o) I 


This completes the proof of the inequality (C.3). 

By combining the inequalities (( .1) and (C.3), we prove the inequality (12). 
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Appendix C.l. Proof of the inequality (C.l) 

We first calculate ||rj^"' ) — (<fa)|| for some q m and q m ~i- Because of the upper 

bound (8) for a small-time evolution, there exists an operator for any P[ q ff~ L 1 ' > such 
that 


||p(9m) _ p(9m-l)/r.\|| - oq m -i/fc (^dt/2) < ~ qm gm A/ ’ ( gm _ 1 ) 

\\ L tm i tm-1 \ UL )W — Z l—n5t/2 11 im-1 11 

< 2 —*— +1 ri 9m ~ } 

— 11 tm-l 11 5 


(C.4) 


where the definition (10) give n5t/2 < 1/2 and in the last inequality, we use the fact 
that the function x^' ?m_l?m_1 ' ) / A: /(l — x) monotonically increases for 0 < x < 1. 

We now define a positive integer 5 q such that 


Qm 2 q m _i T 5q 

for m — 1, 2,..., n. We then obtain 


Qn — 2 n (q l 0 + dq) — d q . 

Because of the condition q n < q, we have to take S q so that it may satisfy the inequality 

a o n a 

2 n (qo +Sq) - 5 q < q, or S q < 2n _ 1 ° - 


Based on this inequality, we choose S q as |_(g — 2 n q 0 )/(2 n — 1)J. By combining the 
inequality (< I) with the definition S q := q m — 2q m _\, we finally obtain 

|i r (9m) _ p(5m-< 2-^ +1 ||r^ m - l) || 

I I’m Cm-1 V •'ll — II C m -1 II 


<4||r<ryll|exp(-log2 


We here notice the equality 
A = 4 exp 


1 q 

c - Qo 


= 4 exp — log 2 


>g /( 2"-l)-g 0 
k 


,g/( 2 n — l) — g 0 
k 


(C.5) 

(C.6) 


because of the dehnitions of £, r 4 and n as in Theorem 3.1 and Eq. (10). We thus prove 
the inequality (C.l) from (C.5) and (C.6). 


Appendix D. Derivation of the inequality (1 ) 


In this section, we consider the norm of 


N 

2=1 

Because A is given by A = a i with 11a*11 = 1 for i e A, we have 

||n^r‘’>ni,|| < ||r<«>|| for \x'-x\<2q, 
||nyr ( « ) n^,|| = o for \x'-x\>2 q 

for any g-local operator T^h We thereby obtain 

|| = \\ U tn^nx+r t )( h i( t )- I ' iq) ) Il trtcc ’ ,r t x'+r t ) 


(D.l) 
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for i e A and VT^ as long as q < r t {\x' — x\ — 1)/2. 

Because hi is supported on the spin i, we can apply Theorem for q 0 = 
hence 


ITT 4 

| [r t x,rtx+r t ) 


< 


,rtx'+rt) || 

inf ||hi(t) -T (9) || 


r^i')£U(q),q<rt.(\x'—x\— 1)/2 


< 8||hj|| \nt] exp 

= 8 \hit] exp 
Therefore, we obtain 


1 

1 ^ 

r t (\x' — x — 1) 

1^1 

L £ \r t 

L 2 J 

)\ 


x' — x\ — 5 
2 £ 


N 


n 


,h n 


< 


{rtXjTtx+rt) 1111 [r t x',r t x'+rt ) 11 — 


Ell" 


A 

[rtx,r t x+rt ) 


i= 1 

< 8N\nt] exp 


\x' — x\ — 5 


rtx 1 ,r t x'+rt) I 


2£ 


which finally yields the inequality (17) with C v = 8 Ne 5 ^ 2 ^\nt] and /i = l/(2£). 


1, and 


(D.2) 


(D.3) 









